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ABSTRACT 

The recent burst in the number of radii measurements of very low-mass stars 
from eclipsing binaries and interferometry of single stars has opened more ques- 
tions about what can be causing the discrepancy between the observed radii and 
the ones predicted by the models. The two main explanations being proposed 
are a correlation between the radius of the stars and their activity levels or their 
metallicities. This paper presents a study of such correlations using all the data 
published to date. The study also investigates correlations between the radii de- 
viation from the models and the masses of the stars. There is no clear correlation 
between activity level and radii for the single stars in the sample. Those single 
stars are slow rotators with typical velocities v ro tsini < 3.0 km A clear corre- 
lation however exists in the case of the faster rotating members of binaries. This 
result is based on the of X-ray emission levels of the stars. There also appears 
to be an increase in the deviation of the radii of single stars from the models as 
a function of metallicity, as previously indicated by Berger et al. (2006). The 
stars in binaries do not seem to follow the same trend. Finally, the Baraffe et 
al. (1998) models reproduce well the radius observations below 0.30-0.35 M , 
where the stars become fully convective, although this result is preliminary since 
almost all the sample stars in that mass range are slow rotators and metallicities 
have not been measured for most of them. The results in this paper indicate that 
stellar activity and metallicity play an important role on the determination of 
the radius of very low-mass stars, at least above O.35M . 

Subject headings: stars: fundamental parameters 
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1. Introduction 

The two most fundamental parameters of a star are its mass and its radius. For a 
given stellar mass, models try to reproduce the radius of the star by implementing the 
best known stellar interior and atmospheric physics. For stars more massive than the Sun, 
convective interior and radiative atmosphere models reproduce well the observations (e.g. 
Andersen 1991, 1997). The atmospheres of these stars are dominated by atomic species and 
their interiors can be closely modeled as an ideal gas. As detailed in the physics of low- 
mass stars review paper by Chabrier & Baraffe (2000), below 1M Q we enter a new physics 
domain where molecular compounds begin to form in the atmospheres of the stars as their 
effective temperature drops, and convection expands to the outer layers of the star, until 
the objects become fully convective below ~ 0.35 M . In addition, the interior of the stars 
becomes denser, and the conditions begin to resemble those of a partially degenerate plasma. 
Therefore the ideal gas equation of state no longer applies. 

Recent measurements of K and M dwarfs radii using both eclipsing binaries and inter- 
ferometry reveal a disagreement between models and observations. As a general trend, the 
measured radii appear to be larger than the predictions by the models by factors as large as 
20-30%, in some cases. Two main hypotheses have been suggested to explain this disagree- 
ment. The first hypothesis suggests that the discrepancies between models and observations 
are caused by differences on the level of activity of the stars (Lopez-Morales & Ribas 2005; 
Ribas 2006). The models by Mullan & McDonald (2001) conclude that the larger radii may 
be caused by the inhibition of convection in stars with strong magnetic fields, or equivalently, 
with high rotational velocities. 

The second hypothesis suggests that the discrepancies between models and observations 
are caused by differences in metallicity (Berger et al. 2006). However, the best current models 
(Baraffe et al. 1998) yield radii values that only differ by about 3% for metallicities between 
-0.5-0.0 dex. This suggests that the models could be missing some important source of 
opacity, as indidated by the authors. A larger range of metallicities in the models is also 
needed. 

The effort of several groups has provided a current sample of 48 radii measurements for 
K and M main sequence stard_|. Berger et al. (2006), Segrensan et al. (2003) and Lane et al. 
(2001) have measured the radii of 14 single nearby dwarfs through interferometry. The radii 



1 RW Lac A (Lacy ct al. 2005) and HS Aur A (Andersen 1991) have not been included in this study, since 
they appear to be slightly evolved (see Torres et al. 2006). The M-dwarf in RXJ2130. 6+4710 (Maxted et 
al. 2004) has been also excluded since, being the companion to an evolved white dwarf, it has most likely 
experienced a common-envelope evolutionary stage. 
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of another 34 stars, members of eclipsing binaries, have been measured by different groups 
(see references in Table 1). Many of the stars in the sample have metallicity estimations 
and measured X-ray fluxes, a good of magnetic activity indicator (Pevtsov et al. 2003). 
Therefore, there are now enough data to begin testing possible correlations between the 
radii of the stars, their activity levels, and their metallicities. 

This paper presents a study of the possible correlations between the deviation of the 
observed stellar radii from the predictions by the models of Baraffe et al. (1998), the X-ray 
luminosity of the stars, and their estimated metallicities. A study of the radii deviation 
from the models as a function of mass is also presented. Section 2 describes the data sample 
used in the analysis. Section 3 shows the analysis of the correlations between the different 
parameters. Finally, the conclusions of the study are summarized in section 4. 

2. Data Sample 

Tables 1 and 2 summarize the main parameters of the current sample of low-mass 
main sequence stars with radii measurements. The tables include the masses and radii of 
the stars, their metallicities, their X-ray-to-bolometric luminosity ratios, their rotational 
velocities, their effective temperatures and their distances. The sample includes 48 objects 
with masses between 0.092 and 0.960 M , covering the entire range of stellar masses between 
1 M and the brown dwarf limit. The precision in the mass and radius measurements varies 
significantly, depending on the apparent brightness of the targets and the observational 
technique used to derive their parameters. In the case of the components of eclipsing binaries 
brighter than V = 14, and some secondaries in F-M and G-M binaries studied by Pont et 
al. (2005) and Bouchy et al. (2005), errors smaller than 2-3% are reached. For all the other 
objects, i.e. fainter binaries and single stars with radii measured through interferometry, the 
error bars are larger, typically 5-10%. In the binaries, the masses and radii of the stars can 
be measured directly from the light and radial velocity curves of the systems. In the case of 
single stars, interferometry provides direct measurements of the stellar radii, but the masses 
need to be derived using external calibrations. Berger et al. (2006), Segransan et al. (2003) 
and Lane et al. (2001) have derived the masses of the stars in their interferometric studies 
using the empirical K-band Mass-Luminosity relation from Delfosse et al. (2000) and the 
Mass-M/f relation by Henry & McCarthy (1993). The authors adopted errors of 5-10% in 
the derived masses to account for photometric and empirical fitting errors. 

Determination of the metallicity of low-mass stars is a difficult task, since as the spectral 
type of the stars increases, i.e. their effective temperature decreases, absorption bands from 
molecules such as TiO, VO, H 2 0, CO and CN begin to appear in the spectra, causing the 
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spectral continuum to be underestimated. Efforts to derive accurate metallicities of M-dwarfs 
are currently underway (see Bonfils et al. 2005, Woolf & Wallerstein 2005 and Bean et al. 
2006). Bonfils et al. (2005) and Woolf & Wallerstein (2005) have measured metallicities for 
most of the single stars in the sample. The values that they derive are summarized in column 
4 of Table 1. Three of the secondary stars in eclipsing binaries, V818 Tau B, OGLE-TR-34 B, 
and OGLE-TR-122 B have also reported metallicities, based on the spectral analysis of their 
more massive primaries (Boesgaard & Friel 1990; Pont et al. 2005). Those values are also 
reported in Table 1. In the case of the M-dwarf eclipsing binaries YY Gem, CU Cnc and GU 
Boo, Torres & Ribas (2002), Ribas (2003) and Lopez-Morales & Ribas (2005) have reported 
abundances close to solar. The metallicity estimates for CU Cnc and GU Boo are based 
on the Galactic dynamics of the systems. However those estimations are not strict enough, 
and have not been considered in the analysis presented in this paper. Torres & Ribas (2002) 
have estimated a metallicity of [Fe/H] = 0.1 ± 0.2 dex for YY Gem, based on metallicity 
measurements of Castor A and B, of which YY Gem is believed to be a companion. Values of 
the metallicity of CM Dra have been derived by several authors using both spectral analysis 
and comparisons with isochrone models (Chabrier & Baraffe 1995; Leggett et al. 2000; Viti 
et al. 2002). They obtain values between [Fe/H]= -0.4 and -1.0 dex. For this study I have 
adopted an average value of [Fe/H] = -0.67 ± 0.20 dex. Finally, metallicity estimations 
of V1061 Cyg, FL Lyr, RW Lac, and HS Aur have been derived by Torres et al. (2006). 
However, those metallicity estimations have not been used, since the are based on model fits 
and the purpose of this work is to compare models to observations. 

Most objects in the sample with apparent magnitudes brighter than V = 14 appear as 
X-ray sources in the ROSAT All-Sky Bright Star Catalog (Voges et al. 1999), with typical 
position offsets between the optical sources and their X-ray counterparts of less than 20-30 
arcsec. The X-ray counts from the ROSAT Catalog have been converted to fluxes using the 
equation 

F x = (5.30HR + 8.31)10- 12 X courits (1) 

derived by Schmitt et al. (1995), where F x is the X-ray flux of each source, HR is its 
hardness ratio and X counts is the number of X-ray counts per second. The units of F x are 
ergs s^ 1 . Those fluxes were then converted to X-ray luminosities using distances from the 
Hipparcos Catalog (ESA 1997) in the case of the single stars and the binaries FL Lyr and 
VI 061 Cyg. The distances to all the other bright binaries were compiled from Chabrier & 
Baraffe (1995), Torres & Ribas (2002), Ribas (2003), Lopez-Morales & Ribas (2005), Lacy et 
al. (2005), and Torres et al. (2006). All the available distances are listed in column 4 of Table 
2. In the case of objects with measured distances, but no X-ray emission detections, X-ray 
luminosity upper limits have been computed assuming a detection limit of 0.005 counts 
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and HR = -1, the minimum value for the hardness ratidj. Those objects are BW3 V38, TrES- 
HerO-07621, 2;mass-j0516288l l+2607387, 2MASS J04463285+1901432, UNSW-TR-2, V1061 
Cyg B, RW Lac and HS Aur. 

Since the stars in the sample cover a wide range of masses, and therefore radii, the X- 
ray-to-bolometric luminosity ratio, Lx/Lboi, gives a more accurate estimate of the relative 
activity level differences between the stars. The bolometric luminosities of the stars have 
been computed from the measured radii and effective temperatures of the objects listed in 
Tables 1 and 2. 

In the case of the binaries, the X-ray luminosities derived above correspond to the 
combined contribution from the two stars. Therefore, it was necessary to device a way 
to estimate the X-ray luminosity of each component. All the binaries in the sample with 
detected X-ray emission have orbital periods of less than 6 days (in fact, only V818 Tau has 
P > 3 days), and orbital solutions consistent with e=0.0, i.e. circular orbits (see references 
in Table 2). The rotational and orbital periods of the stars are therefore expected to be 
synchronized by tidal interaction. Based on this assumption, the rotational velocity of the 
stars can be computed from their radii and the orbital period of the systems. This assumption 
is in fact necessary in the cases where no direct rotational velocity measurements are available 
(see last paragraph in this section). The secondary stars in the binaries are expected to have 
slightly slower rotational velocities than the primaries, and one would therefore intuitively 
expect their magnetic fields to be slightly weaker than the magnetic fields of the primaries. 

X-ray coronal emission has been found to be correlated with stellar rotation in active 
T-Tauri and late- type main sequence single stars (e.g. Bouvier 1990; Fleming et al. 1989). 
However, no correlation between those parameters has been found among binaries (Fleming 
et al. 1989), which are usually faster rotators (v rat sini > 10 kms -1 ). This result, attributed 
to saturation effects, has been corroborated in subsequent studies of L x /L bo i emission of 
G-, K- and M- dwarfs in the Pleiades (Stauffer et al. 1994), the Hyades (Reid et al. 1995; 
Stauffer et al. 1997), and nearby M-dwarfs (Delfosse et al. 1998). These authors find that 
in the cluster M-dwarfs the Lx/L^oi emission saturates at rotational velocities v rot sini > 15 
kms -1 , while for field dwarfs the Lx/L^i emission saturates at v rot sini > 5 kms -1 . In both 
cases the reported saturation level is L x /L bo i ~ 0.003. 

All the components of binaries in the sample with detected X-ray emission have pro- 
jected rotational velocities larger than ~ 7 km s -1 , therefore they are expected to fall within 
the saturated X-ray emission regime, where no correlation between Lx and v rot sini has been 
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found. However, to provide a more complete test of any possible correlation between Lx / t>boi 
and the radii deviation from the models of low-mass stars, I have estimated the X-ray lumi- 
nosity of each binary component considering all the observed correlation cases between Lx 
and v rot sini. These cases are: 

Case 1: No correlation between Lx and v rot sini for all the binaries, as found by (Fleming 
et al. 1989). In this case it has been assumed that each binary component contributes the 
same amount of X-ray flux, i.e. Lxi — Lx2 =Lx/2, where Lx is the total X-ray luminosity 
of the binary measured by ROSAT, and Lxi and Lx2 are, respectively, the luminosities of 
the primary and the secondary. This is a reasonable approximation for the binaries in this 
sample, since most of them have components of similar mass. 

Case 2: A linear correlation between Lx and v rot sini for all the binaries, i.e. Lx/L^i 
oc v rot sini, as found for single main-sequence stars (Fleming et al. 1989). In this case the 
X-ray luminosity of each star in the binaries can be computed using the relations 

Lx = Lxi + Lx2 (2) 

and 

Lxi _ Vrotl ^ 
Lx2 V rot 2 

where v rot i and v rot 2 are their rotational velocities. The resultant expressions for Lxi and 
Lx2 are in this case 



v 



r 



J X1 



l+V r 

and 



Lx (4) 



where v r = v rotl /v rot2 . 

Case 3: A square correlation between L x and v rot sini for all the binaries, i.e. L x /L bo i 
oc (vrotsini) 2 , as found for T-Tauri stars (Bouvier 1990). In this case 
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and the expressions for L X i and L X 2 become 
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and 



^X2 — -t— — o 
1 + 1>; 



The Lx/Lboi obtained in each case are summarized in Table 3. The values of Lx/L^i 
listed in Table 1 correspond to Case 1, since the assumption made in that case is the most 
consistent with the observational results for binaries. 

Finally, Table 2 also includes the estimated rotational velocities for each star. Most of 
the single stars in the sample have rotational velocity estimates, v rot sini, derived directly 
from the rotational broadening of their spectral lines by Delfosse et al. (1998). In this case, 
the reported v rot sini values are upper limits, given by the resolution of the spectra used in 
that work. For the binaries, rotational velocity measurements are only available for some of 
them, i.e. YY Gem A, B (Torres & Ribas 2002), GU Boo A, B (Lopez-Morales & Ribas 2005), 
FL LyrB, V1061 Cyg Ab, B, and RW Lac B (Torres et al. 2006). The rotational velocities 
measurements of these systems are consistent with synchronous rotation in all cases, except 
RW Lac B. RW Lac has an orbital period of 10.3692 days and a reported eccentricity of 
e ~ 0.01 (Lacy et al. 2005), therefore, assuming synchronous rotation mat not be correct. 
The values of v ro tsini gien for the binaries in Table 2 correspond to synchronous rotation 
(v sync sini). They have been calculated from the radii of the stars and the orbital period 
of the binaries published in the literature. In the case of the OGLE binaries reported by 
Bouchy et al. (2005) and Pont et al. (2005), synchronous rotational velocities have been 
computed only for the systems with orbital periods < 5 days (see Note b in Table 2). 



3. Mass, X-ray Luminosity, and Metallicity Correlation with Radius 

The Baraffe et al. (1998) models were able to reproduce, for the first time, the Mass- 
Luminosity and Mass- My relations of low-mass stars. However, further tests including the 
Mass- Radius relation and the effect of metallicity and magnetic activity have been inhibited 
by the lack of observational data. 

The new sample of mass, radius, metallicity and stellar activity data compiled in this 
paper now enables those tests. This section compares the stellar radii predictions of the 
Baraffe et al. (1998) models to the observed radii of the stars as a function of the other 
three parameters, i.e. masses, metallicities and stellar activity. The adopted reference model 
corresponds to an isochrone of age 1 Gigayear, solar metallicity [Fe/H] = 0.0, and mixing 
length a = 1.0 (standard modeljj. 



3 The models yield different radii for M > 0.7 M© when different values of a are used. 1.0 Rq for 1.0 
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Chabrier & Baraffe (1997) and Baraffe et al. (1997) found that variations of a within a 
factor of ~ 2 do not significantly alter the result of their models below ~ 0.6 M & , but those 
variations become important above that mass. More recently, Yildiz et al. (2006) have found 
a very definitive correlation between a and stellar mass for masses greater than 0.77 M & . 
In particular, they find a = 0.99 ± 0.03 for V818 Tau B (~ O.76M ) and larger values of a 
for more massive stars. Based on these results, and the fact that the Baraffe et al. (1998) 
models only consider values of a = 1.0 and 1.9 (this last one only for masses above 0.7 M ), 
the Radius vs. Magnetic Activity and Radius vs. Metallicity correlation tests in §3.2 and 
§3.3 only include the stars in the sample with masses < 0.77 M . 

3.1. Radius versus Mass 

Figure 1 shows the current observational Mass- Radius relation for stars below 1M , 
including all the objects in Table 1. Triangles represent single stars, squares represent 
low-mass secondaries to eclipsing binaries with primaries > 1M , and circles represent the 
components of eclipsing binaries below 1M . The solid line shows the theoretical isochrone 
model from Baraffe et al. (1998), for an age of lGyr, Z = 0.02, and mixing length a = 1.0. 
This is the model used as reference throughout this section. 

The scatter above O.3OM is significant and the errorbars of most of those measurements 
certainly need to be improved. However, there are two clear features that can be emphasized 
at this point. The first one is that the radii measurements for stars below O.3OM seem to 
agree well with the models. O.3O-O.35M is precisely the mass at which stars are believed 
to become fully convective. The second feature is that the situation above O.3OM seems 
completely different. A significant scatter in radius is clear for stars of a given mass. This 
is most evident between 0.35 and 0.70 M , where most measurements clump. The current 
sample includes 27 stars between 0.35-0.70 M , 9 stars above 0.70 M , and 12 stars below 

0.30 Mr,. 



3.2. Radius versus Magnetic Activity 

X-ray emission originating in the magnetically heated corona of the stars is a convenient 
indicator of stellar magnetic activity. As mentioned before, in single stars magnetic activity 



Mq corresponds to a — 1.9 and an age of 4.61 Gyr (Y=0.282). However, the difference in radius with the 
reference model adopted above is only +0.0013 Rq. 
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has been found to be correlated with stellar rotation (Bouvier 1990; Fleming et al. 1989), 
while no correlation apparently exists in the case of binaries (Fleming et al. 1989). Fast 
rotators, such as young stars and tidally spun-up stars in close binaries, generally show 
higher levels of X-ray emission than slow rotating stars of similar mass. 

As shown in Table 2, there are two different populations of rotators in the data sample. 
All the single stars are slow rotators, with v rot sini < 2.9 kms~ l . The components of binaries, 
on the other hand, rotate with velocities ranging from v rot ~ 2.0 to 130 kms^ 1 . Those 
differences in rotational velocities translate into differences in magnetic activity levels of the 
order of a hundred (see column 5 in Table 1). 

The top diagram in Figure 2 shows the fractional deviation of the radii of the stars 
from the lGyr, Z = 0.02 Baraffe et al. (1998) model as a function of Lx/ L bo i for the single 
stars in the sample. The bottom diagram in that figure and the two diagrams in Figure 3 
show the same fractional radii deviation for the stars in binaries. Fig 2-bottom represents 
Case 1 in §2. The two diagrams in Figure 3 represent, respectively, Case 2 (top) and Case 
3 (bottom). Notice that the values of the x-axis in the binary diagrams are approximately 
a hundred times larger than the values in Fig 2-top. This reflects the enhanced magnetic 
astivity levels in the components of binaries. The value of Lx/L bo i measured for GJ551 in 
Fig 2-top is 2.728 10~ 4 ± 6.508 10" 5 , that is, a factor of 20-30 times higher than all the 
other single stars in the sample. That value has been rescaled by a factor of ten to include 
that point in the figure. 

The main conclusions that one can derive from the diagrams in Figures 2 and 3 are that 
there is no apparent correlation between the level of magnetic activity and the deviation of 
the radius from the models in the case of single stars. However, a clear correlation exists 
in the case of the more magnetically active components of binaries. These conclusions are 
endorsed by an statistical analysis of the data samples by computing the significance of the 
correlations in each figure using Pearson's correlation coefficient r. The resultant values of 
r are r = -0.115, for the data in Fig 2-top, and r = 0.890, 0.899, and 0.901, for the data 
in Fig 2-bottom, Fig 3-top, and Fig 3-bottom, respectively. These results indicate that 
the null hypothesis, which states that L x /L bo i and the radii deviation from the models are 
uncorrelated, can be only rejected at the < 30% confidence level for single stars, while for 
the stars in binaries the null hypothesis can be rejected at the > 99% confidence level in all 
the cases. The dotted lines in Figure 2-bottom and Figure 3 show the best least square fits 
to the data, using only the binaries with detected X-ray emission, and taking into account 
the errorbars in both quantities. The results of those least square fits are in each case 



Case 1: (R obs - R mo d)/R m od oc 74.709 L x /L bol ; X 2 = 8.41 
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Case 2: (R obs - R mod )/R mod oc 73.430 L x /L bol ; X 2 = 8.74 
Case 3: (R obs - R mo d)/R m od oc 71.690 L x /L bol ; x 2 = 9.11 

The data used in the least square fits only includes stars with M < 0.77M Q to avoid 
biases introduced by a changing mixing- length parameter (see §3). When the data for stars 
more massive than O.77M (FL Lyr B and V1061 Cyg Ab) are included, the slopes of the 
correlations above change to 65.157 (case 1), 62.511 (case 2) and 59.998 (case 3). The 
Pearson's correlation coefficients become r = 0.833, 0.832, 0.829, and the null hypothesis 
can be still rejected at the same confidence levels obtained above. The dot-dashed lines in 
Figures 2 and 3 show the best least square fits including stars > 0.77 M & . 

3.3. Radius versus Metallicity 

Although there is no clear correlation between the magnetic activity levels and radii 
deviation from the models in the case of single stars (Figure 2-top) the dispersion in radius 
is quite large. That dispersion could be attributed to metallicity differences, as shown below. 

The open circles in Fig. 2-top show stars with [Fe/H] < -0.25, while the filled circles 
show stars with metallicities higher than that value. The crosses correspond to stars with no 
available metallicity measurements. With the exception of GJ380, all stars with [Fe/H] > 
-0.25 show larger radii deviations from the models than stars with [Fe/H] < -0.25. There are 
also signs of a gradient between the metallicity of the stars and their radii. That gradient is 
more clearly illustrated in Figure 4. 

Figure 4 represents the relative deviation of the radii of the stars from the models as a 
function of metallicity. The top diagram corresponds to single stars. The bottom diagram 
shows the metallicity estimations for the low-mass secondaries in binaries with primaries > 
1M , V818 Tau B, OGLE-TR-34 B, and OGLE-TR-122 B, and the M-dwarf binaries YY 
Gem and CM Dra. The radii deviation from the models increases with metallicity in the case 
of single stars, as previously noticed by Berger et al. (2006). However the stars in binaries do 
not appear to follow that same trend. As with the Radius vs. Magnetic Activity correlation 
in §3.2, the conclusions in this section have been statistically checked by computing the 
significant of the correlations using Pearson's correlation coefficient. The values of r are in 
this case r = 0.504 for single stars, and r = -0.209 for the binaries. The null hypothesis 
assuming that [Fe/H] and the relative radius deviation from the models are uncorrelated can 
be rejected at the 90% confidence level in the case of single stars. In the case of the binaries 
the null hypothesis cannot be discarded. The dotted lines in Figure 4 show the best least 
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square fit for each sample, taking into account the errorbars in both quantities. The results 
of those least square fits are 

For Single Stars: (R obs - R m od)/R m od oc 0.20 [Fe/H] ; X 2 = 10.70 
For Binaries: (R obs - R mod ) / R mod oc 0.04 [Fe/H] ; X 2 = 10.18 

In this case, all the stars in both samples have masses below 0.77M Q . 

4. Conclusions 

This paper presents a compilation of all the current mass-radius measurements of main 
sequence stars below 1 M , in addition to all the existing information about their X-ray 
emission levels and metallicities. The goal has been to find any correlation between the 
radii of the stars, their masses, their magnetic activity levels, and their metallicities. The 
stellar radii have been compared to the Baraffe et al. (1998) isochrone model for an age of 
1 Gigayear, solar metallicity (Z=0.02), and mixing length a = 1.0. 

The comparison of the observed Mass- Radius relation to the one predicted by the Baraffe 
et al. (1998) model shows good agreement for stars below O.3O-O.35M . Above those masses, 
the deviation of the measured radii from the values predicted by the models is significant. In 
some cases that difference is as large as 30%. The result for masses below 0.35 M & can be a 
bit misleading, since all the objects in that mass range have rotational velocities v ro tsini < 
10 kms -1 (except OGLE-TR-5B, with v ro tsini = 16.5 kms -1 ). More radii measurements of 
active stars in this mass regime are necessary before we can arrive to any conclusion about 
a radius-magnetic activity correlation for those stars. 

Above O.35M , the sample includes stars with a wide range of magnetic activity levels 
and metallicities. There is in this case a significant scatter in radius for any given stellar 
mass. That scatter can be attributed to differences in the magnetic activity level and the 
metallicity of the stars, as shown in §3. The radii of the stars do not appear to be correlated 
with magnetic activity in the case of slow rotators. Slow rotators are usually single stars, 
old enough to have had time to slow down rotationally. The typical rotational velocities of 
these objects are less than 5 km s -1 . However, there is a clear correlation between radius 
and magnetic activity in the case of faster rotating, more magnetically active, stars. It has 
been postulated in previous papers (e.g. Lopez-Morales & Ribas 2005) that the significant 
areal coverage of active regions (spots) on the surface of magnetically active stars changes 
the overall photospheric temperature of the stars; an effect that they could compensate 
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by increasing their radii to conserve their total radiative fluxes. Torres et al. (2006) have 
compared the radii of two fast rotators (V1061 Cyg Ab and FL Lyr B) and two slow rotators 
(RW Lac A and HS Aur A) of similar mass and find the fast rotators to have larger radii than 
what the models predict. The theoretical work by Mullan & McDonald (2001) concludes 
that inhibition of convection in fast rotating stars may be the cause of their larger radii. 
This might be the case not only for low-mass stars in close binaries, but also for young single 
stars and T Tauri objects, which are usually very active. 

Metallicity also appears to be playing a role in the radii of the stars. There is a clear 
correlation between the size of the stars and how metal rich they are. This correlation, 
already noted by Berger et al. (2006), appears to apply to single stars, but not to the 
components of binaries. However, this last conclusion is still premature, since only a handful 
of the low-mass stars in binaries have metallicity estimations. Furthermore, most of those 
metallicity estimations are indirect. Direct metallicity measurements of the components of 
low-mass binaries are needed to determine any correlation with the radii of the stars. An 
increasing abundance of metals would have the effect of enhancing the number density of 
molecular compounds in the atmospheres of the stars, making it harder for the radiation to 
escape. The response of the stars might be similar to the case of high magnetic activity, i.e. 
they could increase their radii to conserve their radiative fluxes. 
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Table 1: Mass, radius, metallicity and X-ray-to-Bolometric luminosity ratio of the stars in 
the sample 



Star 


M 


R 


[Fe/H] 


Lx/L boi 


Source 1 




(Mq) 


(*©) 


(dcx) 







GJ 15 A 0.4040±0.0404 0.379±0.006 -0.46±0.20 

GJ 514 0.5260±0.0526 0.611±0.043 -0.24±0.20 

GJ 526 0.5020±0.0502 0.493±0.033 -0.31±0.20 

GJ 687 0.4010±0.0401 0.492±0.038 0.11±0.20 

GJ 752 A 0.4840±0.0484 0.526±0.032 -0.05±0.20 

GJ 880 0.5860±0.0586 0.689±0.044 -0.04±0.20 

GJ 205 0.631±0.031 0.702±0.063 0.21±0.13 

GJ 887 0.503±0.025 0.491±0.014 -0.22±0.09 

GJ 191 0.281±0.014 0.291±0.025 -0.90±0.20 

GJ 551 0.123±0.006 0.145±0.011 

GJ 699 0.158±0.008 0.196±0.008 -0.50±0.30 

GJ 411 0.403±0.020 0.393±0.008 -0.42±0.07 

GJ 380 0.670±0.033 0.605±0.020 -0.03±0.14 

GJ 105 A 0.790±0.039 0.708±0.050 

CM Dra A 0.2307±0.0010 0.2516±0.0020 -0.67±0.20 

CM Dra B 0.2136±0.0010 0.2347±0.0019 -0.67±0.20 

YY Gem A 0.5992±0.0047 0.6191±0.0057 0.10±0.20 

YY Gem B 0.5992±0.0047 0.6191±0.0057 0.10±0.20 

CU Cnc A 0.4333±0.0017 0.4317±0.0052 

CU Cnc B 0.3980±0.0014 0.3908±0.0094 

GU Boo A 0.610±0.007 0.623±0.016 

GU Boo B 0.599±0.006 0.620±0.020 

BW3 V38 A 0.44±0.07 0.51±0.04 

BW3 V38 B 0.41±0.09 0.44±0.06 

TrES-HerO-07621 A 0.493±0.003 0.453±0.060 

TrES-HerO-07621 B 0.489±0.003 0.452±0.050 

2MASS J05162881+2607387 A .... 0.787±0.012 0.788±0.015 

2MASS J05162881+2607387 B .... 0.770±0.009 0.817±0.010 

2MASS J04463285+1901432 A 0.47±0.05 0.57±0.02 

2MASS J04463285+1901432 B 0.19±0.02 0.21±0.01 

UNSW-TR-2 A 0.529±0.035 0.641±0.05 

UNSW-TR-2 B 0.512±0.035 0.608±0.06 

V818 Tau B 0.7605±0.0062 0.768±0.010 0.13±0.02 

FL Lyr B 0.960±0.012 0.962±0.028 

V1061 Cyg Ab 0.9315±0.0068 0.974±0.020 

V1061 Cyg B 0.925±0.036 0.870±0.087 

RW Lac B 0.870±0.004 0.964±0.004 

HS Aur B 0.879±0.017 0.873±0.024 

OGLE-TR-5 B 0.271±0.035 0.263±0.012 

OGLE-TR-6 B 0.359±0.025 0.393±0.018 

OGLE-TR-7 B 0.281±0.029 0.282±0.013 

OGLE-TR-18 B 0.387±0.049 0.39±0.04 

OGLE-TR-34 B 0.509±0.038 0.435±0.033 0.32±0.31 

OGLE-TR-78 B 0.243±0.015 0.24±0.013 

OGLE-TR-106 B 0.116±0.021 0.181±0.013 

OGLE-TR-120 B 0.47±0.04 0.42±0.02 

OGLE-TR-122 B 0.092±0.009 0.120±0.019 0.15±0.36 

OGLE-TR-125 B 0.209±0.033 0.211±0.027 



2.029E-05±3. 
1.509E-05±6. 
2.653E-06±4. 
9.433E-06±2. 
4.088E-06±7. 
6.296E-06±2. 
1.500E-05±4. 
5.245E-06±1. 
9.904E-06±2. 
2.728E-04±6. 
2.880E-06±3. 
7.169E-06±1. 
8.954E-06±2. 
2.082E-06±8. 



627E-06 
070E-06 
583E-07 
027E-06 
631E-07 
627E-06 
533E-06 
762E-06 
427E-06 
508E-05 
336E-07 
384E-06 
063E-06 
738E-07 



1,15 
1,15 
1,15 
1,15 
1,15 
1,15 
2,15 
2,15 
2,15 
2 

3,17 
3,15 
3,15 
3 



6.518E-04±1.037E-04 
7.858E-04±1.259E-04 
1.326E-03±1.460E-04 
1.326E-03±1.460E-04 
8.457E-04±1.755E-04 
1.079E-03±2.305E-04 
1.191E-03±4.991E-04 
1.347E-03±5.688E-04 

< 1.08E-03 

< 1.53E-03 

< 1.90E-04 

< 1.35E-04 

< 7.70E-04 

< 7.48E-04 

< 1.95E-03 

< 2.44E-02 



4,16 

4,16 

5 

5 

6 

6 

7 

7 



9 

9 

10 

10 

11 

11 

12 

12 



4.025E-04±7.738E-05 
2.819E-04±5.975E-05 
2.373E-04±7.190E-05 

< 1.72E-05 

< 1.07E-05 

< 4.66E-06 



5,18 
19,20 
19 
19 

19,21 

19,22 

14 

14 

14 

14 

14,15 

13 

13 

13 

13 

13 



Bcrger et al. (2006), [2] Segransan et al. (2003), [3] Lane et al. (2001), [4] Metcalfe et al. (1996), 
[5] Torres & Ribas (2002), [6] Ribas (2003), [7] Lopez-Morales & Ribas (2005), [8] Maceroni & Montalban 
(2004), [9] Creevcy et al. (2005), [10] Bayless & Orosz (2006), [11] Hebb et al. (2006), [12] Young et al. 
(2006), [13] Pont et al. (2005), [14] Bouchy et al. (2005), [15] Bonfils et al. (2005), [16] average of values 
derived by Chabrier & Baraffe (1995), Legget et al. (2000) and Viti et al. (2002), [17] Jones et al. (2002), 
[18] Boesgaard & Friel (1990), [19] Torres et al. (2006), [20] Popper et al. (1986), [21] Lacy et al. (2005), 
[22] Andersen (1991). 



Table 2: Rotational velocity, effective temperature and distance of the stars in the sample 



Star 


v ro tsini 


T eff 


Distance 


Source 3, 






(K) 


(parsecs) 





1.4E-05 
9.7E-05 
4.7E-05 
1.7E-05 
6.5E-05 
7.1E-05 
4.7E-05 
8.2E-06 
1.13E-05 
9.8E-06 
8.3E-06 
5.4E-06 
1.6E-05 
5.6E-05 



GJ 15 A < 2.9 3747 ± 100 3.57 ± 

GJ 514 < 2.9 3377 ± 100 7.63 ± 

GJ 526 < 2.9 3662 ± 100 5.43 ± 

GJ 687 < 2.8 3142 ± 100 4.53 ± 

GJ 752 A < 2.6 3390 ± 120 5.87 ± 

GJ 880 < 2.8 3373 ± 100 6.89 ± 

GJ 205 < 2.9 3520 ± 170 5.69 ± 

GJ 887 ■ ■ ■ 3626 ± 56 3.29 ± 

GJ 191 • • ■ 3570 ± 156 3.92 ± 

GJ 551 ••■ 3042 ± 117 1.29 ± 

GJ 699 < 2.8 3163 ± 65 1.82 ± 

GJ 411 < 2.9 3570 ± 42 2.55 ± 

GJ 380 2.8 3950 ± 161 4.87 ± 

GJ 105 A 4714 ± 67 7.21 ± 

CM Dra A 10.03 ± 0.08 3360 ± 100 15.9 ± 0.1 

CM Dra B 9.35 ± 0.08 3320 ± 100 15.9 ± 0.1 

YY Gem A 38.40 ± 0.40 3820 ± 100 14.90 ± 1.0E-04 

YY Gem B 38.40 ± 0.40 3820 ± 100 14.90 ± 1.0e-04 

CU Cnc A 7.86 ± 0.09 3160 ± 150 12.81 ± 5.0E-03 

CU Cnc B 7.11 ± 0.17 3125 ± 150 12.81 ± 5.0E-03 

GU Boo A 64.37± 1.65 3920 ± 130 140 ± 8 

GU Boo B 64.06± 2.07 3810 ± 130 140 ± 8 

BW3 V38 A 129.5 ± 10.2 3500 ~ 400 

BW3 V38 B 111.7 ± 15.3 3448 ~ 400 

TrES-HerO-07621 A 20.28 ± 2.67 3500 118 ± 13 

TrES-HerO-07621 B 20.24 ± 2.24 3395 118 ± 13 

2MASS J05162881+2607387 A .... 15.31 ± 0.29 4200 ± 200 753 ± 34 

2MASS J05162881+2607387 B .... 15.87 ± 0.19 4154 ± 200 753 ± 34 

2MASS J04463285+1901432 A 46.02 ± 1.61 3320 ± 150 ~ 542 

2MASS J04463285+1901432 B 16.95 ± 0.81 2910 ± 150 ~ 542 

UNSW-TR-2 A 15.21 ± 1.19 • ■ ■ 169 ± 14 

UNSW-TR-2 B 14.42 ± 1.42 • ■ ■ 169 ± 14 

V818 Tau B 6.90 ± 0.09 4220 ± 150 46.8 ± 0.4 

FL Lyr B 22.28 ± 0.65 5300 ± 100 130.04 ± 0.025 

V1061 Cyg Ab 20.97 ± 0.43 5300 ± 150 160.00 ± 0.03 

V1061 Cyg B ■ ■ • 5670± 150 160.00 ± 0.03 

RW Lac B b ■ ■ • 5560 ± 150 190 ± 10 

HS Aur B ■ ■ • 5200 ± 75 99.5 ± 0.05 

OGLE-TR-5 B 16.39 ± 0.75 

OGLE-TR-6 B 4.37 ± 0.20 

OGLE-TR-7 B 5.24 ± 0.24 

OGLE-TR-18 B 8.77 ± 0.90 

OGLE-TR-34 B b 

OGLE-TR-78 B b 

OGLE-TR-106 B 3.59 ± 0.26 

OGLE-TR-120 B b 

OGLE-TR-122 B b 

OGLE-TR-125 B 2.00 ± 0.26 



1,2,3 

1,2,3 

1,2,3 

1,2,3 

1,2,3 

1,2,3 

1,4,3 

4,3 

4,3 

4,3 

1,4,3 

1,4,3 

1,21,3 

21,3 



5,6 
5,6 
7 
7 



9 

9 

10 

10 

11 

11 

12 

12 

13 

13 

14 

14 



7,3 
15,16 
15 
15 

15,19 

15,20,3 

17 

17 

17 

17 



18 



18 



Q [l] Delfosse et al. (1998), [2] Berger et al. (2006), [3] Hipparcos Catalog (ESA 1997), [4] Segransan et al. 
(2003), [5] Metcalfe et al. (1996), [6] Chabrier & Baraffe (1995), [7] Torres & Ribas (2002), [8] Ribas (2003), 
[9] Lopez-Morales & Ribas (2005), [10] Maccroni & Montalban, [11] Creevey et al. (2005), [12] Bayless & 
Orosz (2006), [13] Hebb et al. (2006), [14] Young et al. (2006), [15] Torres et al. (2006), [16] Popper et al. 
(1986), [17] Bouchy et al. (2005), [18] Pont et al. (2005), [19] Lacy et al. (2005), [20] Andersen (1991), [21] 
Ramirez & Melendez (2004) 

b ow-mass secondaries in binaries with orbital periods > 5 days. Three of these binaries show eccentricities 
c > 0.0 (Pont et al. 2005). Therefore, the rotational velocities of the stars can not be computed from the 
assumption that the systems are tidally locked. This is consistent with the theoretical result from Zahn 
(1989), that predicts larger orbital circularization timescales for binaries with small mass ratios. Pont et 
al. (2005) conclude from their data that the typical circularization period for these systems is 5 days. The 
binary OGLE-TR-125A-B has a period of 5.3039 days, but eccentricity e = 0.0 ± 0.01. therefore its orbit is 
circularized. 



Table 3: X-ray-to-bolometric luminosity ratio of the stars in each binary with detected X-ray 
emission. Case 1 corresponds to no correlation between Lx and w r0 (Sini, Case 2 corresponds 
to L x oc v rot sim, and Case 3 corresponds to L x oc (v rot smi) 2 . 



Case 1 Case 2 Case 3 

Star Lx/L bo i Lx/L boi Lx/L bol 

CM Dra A 6.518E-04 ± 1.037E-04 6.745E-04 ± 1.073E-04 6.971E-04 ± 1.111E-04 

CM Dra B 7.858E-04 ± 1.259E-04 7.585E-04 ± 1.216E-04 7.313E-04 ± 1.175E-04 

YY Gem A 1.326E-03 ± 1.460E-04 1.326E-03 ± 1.462E-04 1.326E-03 ± 1.470E-04 

YY Gem B 1.326E-03 ± 1.460E-04 1.326E-03 ± 1.462E-04 1.326E-03 ± 1.470E-04 

CU Cnc A 8.457E-04 ± 1.755E-04 8.877E-04 ± 1.846E-04 9.295E-04 ± 1.943E-04 

CU Cnc B 1.079E-03 ± 2.305E-04 1.025E-03 ± 2.195E-04 9.719E-04 ± 2.096E-04 

GU Boo A 1.191E-03 ± 4.991E-04 1.194E-03 ± 5.009E-04 1.197E-03 ± 5.039E-04 

GU Boo B 1.347E-03 ± 5.688E-04 1.344E-03 ± 5.681E-04 1.341E-03 ± 5.688E-04 

V818 Tau B 4.025E-04 ± 7.738E-05 3.706E-04 ± 7.140E-05 3.392E-04 ± 6.578E-05 

FL Lyr B 2.819E-04 ± 5.975E-05 2.417E-04 ± 5.148E-05 2.031E-04 ± 4.408E-05 

V1061 Cyg Ab 2.373E-04 ± 7.190E-05 1.786E-04 ± 5.416E-05 1.266E-04 ± 3.859E-05 
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Mass (M Q ) 



Fig. 1. — Current observational Mass- Radius relation for stars below 1M Q (objects in Table 
1). The top figure shows all the data from low-mass secondaries to eclipsing binaries with 
primaries > 1M© (squares), and the components of eclipsing binaries below 1M (circles). 
The bottom figure shows all the measurements from single stars (triangles). The solid line 
in both figures represents the theoretical isochrone model from Baraffe et al. (1998), for an 
age of lGyr, Z = 0.02, and mixing length a = 1.0 (standard model). 
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Fig. 2. — Top - Fractional deviation of the radius of the single stars in the sample from the 
lGyr, Z = 0.02, a = 1.0 Baraffe et al. (1998) model as a function of Lx/ L^. Filled circles 
show stars with [Fe/H] > -0.25 dex and open circles show stars with [Fe/H] < -0.25 dex. 
Stars with no available metallicity estimations are shown as crosses. Bottom - Fractional 
deviation of the radius of the stars in binaries from the lGyr, Z = 0.02, a = 1.0 Baraffe et 
al. (1998) model as a function of L x /L ho i for Case 1, i.e. assuming no correlation between 
Lx/Lboi and iwsini. The open circles show the component of binaries with X-ray emission 
detected by ROSAT, the crosses indicate L x /L bo i upper limits for the component of binaries 
with no X-ray emission detected. The dotted line shows the best least square fit to the data 
(open circles, M < O.77M ). The dot-dashed line shows the best least square fit to the data 
including stars with M > O.77M . The dashed line in both figures indicates the model's 
zero-deviation baseline. 
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Fig. 3. — Top - Fractional deviation of the radii of the binary stars in the sample with 
detected X-ray emission from the lGyr, Z = 0.02, a = 1.0 Baraffe et al. (1998) model 
as a function of Lx/Lb Q i for Case 2, i.e. assuming Lj/Lw oc t> rot sini. Bottom - Same as 
top figure for Case 3, i.e. assuming Lx/Lboi oc (v rot sini) 2 . The dotted lines show the best 
least square fits to the data with M < O.77M . The dot-dashed lines show the best least 
square fit to the data including stars with M > 0.77M Q . The dashed lines show the model's 
zero-deviation baseline. 
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Fig. 4. — Fractional deviation of the radius of stars in the sample from the lGyr, Z = 0.02, 
a = 1.0 Baraffe et al. (1998) model as a function of metallicity. The Top diagram shows 
the single stars. The Bottom diagram shows the binary low-mass secondaries V818 Tau B, 
OGLE-TR-34 B, and OGLE-TR-122 B, and the M-dwarf binaries YY Gem and CM Dra. 
The dotted lines show the best least square fits to the data. The dashed lines show the 
model's zero-deviation baseline. 



